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ABSTRACT 



The strongest evidence for dark energy comes presently from geometric tech- 
niques such as the supernova distance-redshift relation. By combining the mea- 
sured expansion history with the Friedmann equation one determines the energy 



density and its time evolution, hence the equation of state of dark energy. Be- 
—i- . cause these methods rely on the Friedmann equation which has not been inde- 

pendently tested it is desirable to find alternative methods that work for both 
O . general relativity and other theories of gravity. 

Assuming that sufficiently large patches of a perturbed Robertson- Walker 



spacetime evolve like separate Robertson- Walker universes, that shear stress is 
unimportant on large scales and that energy and momentum are locally con- 
served, we derive several relations between long- wavelength metric and matter 
perturbations. These relations include generalizations of the initial-value con- 



straints of general relativity. For a class of theories including general relativity 
we reduce the long-wavelength metric, density, and velocity potential perturba- 
tions to quadratures including curvature perturbations, entropy perturbations, 
and the effects of nonzero background curvature. When combined with the ex- 
pansion history measured geometrically, the long-wavelength solution provide a 
test that may distinguish modified gravity from other explanations of dark energy. 

Subject headings: cosmology:theory,gravitation 



1. Introduction 

Current evidence for dark energy is based on two key assumptions: 
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1. The Cosmological Principle holds, i.e. on large scales the matter distribution and its 
expansion are homogeneous and isotropic and the spacetime geometry is Robertson- 
Walker; and 

2. The cosmic expansion scale factor a(t) obeys the Friedmann equation, which may be 

written 

2 _/lda\ 2 22 8irGa 2 
H = (^J =aH =—^-p(a)-K, (1) 

where K is the spatial curvature with units of inverse length squared (we set c = 1) and r 
is conformal time, related to cosmic proper time t by dt = o(r)dr. Usually the Friedmann 
equation is applied indirectly through an expression for the angular-diameter or luminosity 
distance; these formulae depend crucially on the expression for r(a) obtained by integrating 
equation (1). 

The Cosmological Principle is more general than general relativity (GR) and it is 
amenable to direct observational test through measurements of distant objects such as galax- 
ies, Type la supernovae, and the microwave background radiation. The Friedmann equation, 
however, is equivalent to one of the Einstein field equations of GR applied to the Robertson- 
Walker metric and, so far at least, it has not been independently tested. Instead, the 
Friedmann equation is used by astronomers in effect to deduce p(a) including dark energy 
through measurements of the Hubble expansion rate H(a). 

If the evidence for dark energy is secure, there are four possible explanations: 

1. The dark energy is a cosmological constant or, equivalently, the energy density and 
negative pressure of the vacuum (Gliner 1966; Zel'dovich 1967). 

2. The dark energy is some other source of stress-energy, for example a scalar field with 
large negative pressure (Ratra & Peebles 1988; Steinhardt et al. 1999). 

3. General relativity needs to be modified and we cannot use it to deduce the existence 
of either a cosmological constant or an exotic form of energy (Dvali et al. 2000; Lue et 
al. 2004). 

4. General relativity is correct but our understanding of it is not; e.g. long-wavelength 
density perturbations modify the Friedmann equation (Kolb et al. 2005a,b). 

It would be interesting to find observational tests that can distinguish among these 
possibilities. The first case can be tested by measuring the equation of state parameter 
w = p/p on large scales; a cosmological constant has w — — 1. This test can be made using 
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methods such as the supernova distance-redshift relation (Riess et al. 1998; Perlmutter et al. 
1999) and baryon acoustic oscillations (Eisenstein et al. 2005) whose interpretation relies on 
the Friedmann equation. Such methods are called geometric (they measure the large-scale 
geometry of spacetime) or kinematic (they rely on energy conservation and on the first time 
derivative of the expansion scale factor). 

It is more difficult to test the third possibility, namely that dark energy represents a 
modification of GR rather than (or in addition to) a new form of mass-energy Redshift- 
distance tests invoke the Friedmann equation of GR or its equivalent in other theories in 
order to determine the abundance and equation of state of dark energy. Such tests do not 
work without a formula relating a(r) to the cosmic energy density and pressure. While 
specific alternative models can be tested, it would be desirable to have general cosmological 
tests independent of the Friedmann equation. 

The evolution of density perturbations has been proposed as an independent test of 
dark energy and general relativity (Linder 2005; Ishak et al. 2005). For example, weak 
gravitational lensing measurements, the evolution of galaxy clustering on large scales, and 
the abundance of rich galaxy clusters all have some sensitivity to the gravitational effects of 
dark energy at redshifts z < 1. These methods may be called dynamic because the evolution 
equations for the perturbations are at least second order in time. 

This paper examines what the linear growth of metric, density, and velocity pertur- 
bations can tell us about gravity. By generalizing previous work on evolution of separate 
universes (Wands et al. 2000; Gordon 2005 and references therein), we will show that in GR, 
on length scales larger than the Jeans length (more generally, on scales large enough so that 
spatial gradients may be neglected in the equations of motion), the evolution of the metric 
and density perturbations of a background Robertson- Walker universe can be determined 
from the Friedmann equation and local energy-momentum conservation. Generalizing this 
to arbitrary theories of gravity, under certain conditions explained in this paper, the evolu- 
tion of Robert son- Walker spacetimes (assuming they are solutions of the gravitational field 
equations) combined with local energy-momentum conservation is sufficient to determine the 
evolution of long-wavelength perturbations of the metric and matter variables. 

Because GR is so fully integrated into most treatments of cosmological perturbation 
theory, and we wish to test gravitation more generally, it is worth recalling its basic elements: 

1. Spacetime is describable as a classical 4-dimensional manifold with a metric locally 
equivalent to Minkowski. This is generalizable to higher dimensions where matter 
fields reside on a three-dimensional spatial brane. 

2. Special relativity holds locally. In particular, energy-momentum is locally conserved. 
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This is more general than GR. 

3. The weak equivalence principle holds, i.e. freely-falling bodies follow spacetime geodesies. 
This is more general then GR. 

4. The metric is the solution to the Einstein field equations subject to appropriate initial 
and boundary conditions. This is uniquely true in GR. 

The first three ingredients are assumed by most viable theories of gravitation applied on 
cosmological scales. They are assumed to be correct throughout this paper. The Einstein 
field equations are not assumed to hold except where explicitly stated below. We assume 
throughout that the universe is approximately (or, for some calculations, exactly) Robertson- 
Walker. 



2. Robertson- Walker Spacetimes and their Perturbations 

The general Robertson- Walker spacetime is specified by giving a spatial curvature con- 
stant K (with units of inverse length squared) and a dimensionless scale factor a(r) (nor- 
malized so that a = 1 today). Let us assume that the evolution of a(r) depends on K and 
on the properties of the matter and energy filling the universe. Applying the Cosmological 
Principle, the matter must behave as a perfect fluid at rest in the comoving frame. The 
pressure of a perfect fluid may be written p(p, S) where p is the proper energy density and 
S is the comoving entropy density in the fluid rest frame. The line element may thus be 
written 

ds 2 = cl 2 {t, K, S) [-dr 2 + d X 2 + r 2 { X , K)dtt 2 ] (2) 

where r(x,K) = K~ 1 / 2 sm.(K 1 / 2 x) for K > (and is analytically continued for K < 0) 
and dfi 2 = d6 2 + sin 2 6 d(p 2 . Initially we make no assumption about the dynamics except 
that the evolution of the scale factor depends only on the geometry (K) and composition 
(S) of the (3+l)-dimensional universe and not, for example, on parameters describing extra 
dimensions. 1 Only later will the Friedmann equation be assumed to determine the exact 
form of a(r, K, S). 



x It would be straightforward to add such parameters to the argument list of a(r, K, S) and then perturb 
them. 
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2.1. Curvature and Coordinate Perturbations 

Metric perturbations are obtained by comparing two slightly different spacetimes. We 
consider two homogeneous and isotropic Robert son- Walker spacetimes differing only by their 
(spatially homogeneous) spatial curvature. The first spacetime has spatial curvature K; the 
second one has spatial curvature K{\ + 8k) where 8k is a small constant. We write the 
metric of the second spacetime as a perturbation of the first, as follows. First, the angular 
radius may be Taylor-expanded to first order in 8k to give 

r 2 ( X , K + K8 K ) = r 2 ( X , K)(l - 8 K ) + X r( X , 4K)5 K , (3) 

while a(r, K + K8k, S) = a(r, K, S) + <5x(<9a/<91n K). To simplify the appearance of the line 
element we change variables r — ► r + a(r) and X ~^ x(l ~~ K ) where a and k are assumed to 
be first order in 8k- In these new coordinates, to first order in 8k the second spacetime has 
line element 

ds 2 = a 2 (T,K,S)(l + 2na + 25 K ^^] (4) 



x 



{-(1 + 2a)dr 2 + (1 - 2K)d X 2 + [(1 - 5 K )r 2 ( X , K) + (8 K - 2K) X r( X , *K)] ^} , 



where Ti(r,K,S) = din a/ dr. This line element describes a perfectly homogeneous and 
isotropic Robertson- Walker spacetime with spatial curvature K(l + 8k)- However, for ap- 
propriate choices of a and k (i.e., the appropriate coordinate transformation), it takes pre- 
cisely the same form as a perturbed Robert son- Walker spacetime with background spatial 
curvature K and perturbation ^, 

ds 2 = a 2 {r, K,S){-(1 + 2^)dr 2 + (1 - 2^) [d X 2 + r 2 { X , K)dtt 2 } } , (5) 

provided that the following three conditions hold: 

k=^8 k , 2V = a + K, ^ = K ( 1 - 2 §]^)-' Ha - ( 6 ) 

In other words, a perturbed Robertson- Walker spacetime whose metric perturbation \P(t) 
is spatially homogeneous is identical to an unperturbed Robertson- Walker spacetime repre- 
sented by a perturbed coordinate system. 2 

This equivalence is significant because it suggests that long-wavelength curvature per- 
turbations (for which \1/ is effectively independent of spatial position) should evolve like 



2 In §4 we will generalize equation (5) to the case of two distinct potentials for the time- and space- parts 
of the metric. 
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patches of a Robert son- Walker spacetime (Wands et al. 2000; Gordon 2005 and references 
therein), whose dynamics is more general than general relativity. So far we have assumed 
that the Cosmological Principle holds but we have not assumed the validity of the Einstein 
field equations of general relativity. However, if we know how a perfect Robertson- Walker 
spacetime evolves, the evolution of long-wavelength curvature perturbations follows. 

To make further progress let us assume the form of a(r, K, S), which requires specifying 
a theory of gravity. The simplest choice is to assume the validity of the Friedmann equation, 
which gives 

"8vr 



r(a, K, S) 

Using this, one finds 

d In a 



1-1/2 



d\nK 



—Ga 4: p(a, S) — Ka 2 



(dr/d\nK) a , s __1 
2 



da . 



dr' 



r>3 (dr/d\na) K , s 2 J H 2 (r',K,S) 

We also require conservation of energy, which may be written 

d 



d In a 



p(a,S) = -3[p + p(p,S)] . 



When combined with the Friedmann equation, this implies 



(7) 



(9) 



7 = 4%Ga 2 (p + p) = H 2 + K -H 



-(l + w)(H 2 + K) 



(10) 



where H = &H{r, K, S)/dr. Combining equations (6), (8), and (10) yields a relation between 
the long-wavelength potential \P and the curvature perturbation 3 k: 

H 2 d /a 2 *' 



K 



7a 2 dr \ H 



;n) 



For long wavelengths k may depend on the wavenumber k but cannot depend on r, so 
8k/ dr = 0, which implying 

' J/ °'"^\ {jj + 3 (x + c 2 w )H^ + 3(c 2 w -w)H 2 ^- (2 + 3w + 3c 2 w )K^ = 0, (12) 



7 d 



7a 2 dr \ H 



where 



w 



p(p, S) 



dp\ 
d P/s 



(13) 



For spatially uniform, iscntropic perturbations with vanishing shear stress, k reduces to the £ variable 
of Bardcen ct al. (1983). In other cases the two variables differ, with k being simpler in both its dynamics 
and interpretation. 
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By comparing two different Friedmann-Robertson- Walker models we have arrived at a 
second-order differential equation for the metre perturbation \l/. The only dynamical equa- 
tions assumed have been the Friedmann equation and energy conservation in a homogeneous 
and isotropic universe. We have not made use of the perturbed Einstein or fluid equations. 
Nevertheless, as we show next, equation (12) is identical with the dynamical evolution equa- 
tion for long- wavelength curvature perturbations obtained using the perturbed fluid and 
Einstein equations of general relativity. 



2.2. Linear Cosmological Perturbations in General Relativity 

Cosmological perturbation theory has been well studied (e.g. Lifshitz 1946; Bardeen 
1980; Kodama & Sasaki 1984; Hwang & Noh 2002). Nonetheless, the curvature and en- 
tropy variables relevant for long-wavelength density perturbations differ from those given 
previously in the literature, so a brief summary is presented here. 

In the conformal Newtonian gauge (Mukhanov et al. 1992; Ma & Bertschinger 1995), 
the metric of a perturbed Robertson- Walker spacetime with scalar perturbations may be 
written as a generalization of equation (5), 

ds 2 = a 2 (T){-(l + 2^)dr 2 + (l-2^)[d X 2 + r 2 (x,K)dn 2 }} . (14) 

Here, §(x l ,r) and ^(x l ,r) are small-amplitude gravitational potentials. The dependence of 
the scale factor on K and S is suppressed because we are considering now a single universe 
with unique values of these parameters and are introducing perturbations only through the 
potentials. 

For scalar perturbations the stress-energy tensor components may be written in terms 
of spatial scalar fields 5p, u, and n, as follows, 

T° = -ip + Sp), (15a) 

T°, = -(p + p)V lU , (15b) 

^(p + p)(Vv,-i ( 



T\ = S i j (p + Sp) + -(p + p)(v i V j --8 i j A)7r : (15c) 



where Vj and V* are the three-dimensional covariant derivative for the spatial line element 
dx 2 + r 2 dQ 2 while A = V*Vj. Unless stated otherwise, all variables refer to the total stress- 
energy summed over all components. The unperturbed energy density and pressure are p(r) 
and p(t), respectively, while 5p and 5p are the corresponding perturbations measured in the 
coordinate frame. 



In the scalar mode, all perturbations to the metric and the stress-energy tensor arise 
from spatial scalar fields and their spatial gradients. For example, the energy flux is a 
potential field with velocity potential u(x 1 ,t). Similarly, the shear stress follows from a 
shear stress (viscosity) potential 7i(x 1 ,t) (defined as in Bashinsky & Seljak 2004). For an 
ideal gas, n = 0. Equations (15) are completely general for the scalar mode. Stress-energy 
perturbations that arise from divergenceless (transverse) vectors contribute only to the vector 
mode, while divergenceless, trace-free tensors contribute only to the tensor mode. Vector 
and tensor modes are ignored in this paper. 

Combining the perturbed Einstein field equations (Bertschinger 1996) with equations 
(10) and (13) gives the following second-order partial differential equation for the linear 
evolution of the potential \I/: 



7 d 



H 2 d fa 



■fa 2 dr \H 



—^ 



c 2 w A^ = 7 



Sp - c 2 w 5p 3 d 



p + p Hdr 



+ ^(^-)+ A - 



(16) 



Equation (16) is exact in linear perturbation theory of general relativity and is fully general 
for scalar perturbations. It was given in another form by Hwang & Noh (2002). Aside from 
the right-hand side and the sound speed term c^A\I/, it agrees exactly with equation (12), 
whose derivation did not assume the validity of the perturbed Einstein field equations. The 
right-hand side of equation (16) includes entropy perturbations proportional to Sp — c^Sp 
and shear stress perturbations, both of which were absent in the Robertson- Walker models 
of the previous subsection. 

In a perfectly homogeneous and isotropic universe the shear stress must vanish, and a 
correct treatment of shear stress requires going beyond the Friedmann and energy conserva- 
tion equations. In the universe at low redshift, the shear stress due to photons, neutrinos, 
and gravitationally bound structures is orders of magnitude smaller than the mass density 
perturbations. Unless the dark energy is a peculiar substance with large shear stress, we 
may neglect 7r in the equation of motion for \l/. 

The sound wave term represents the effect of pressure forces in resisting gravitational 
instability. For a perturbation of comoving wavenumber k, — c^A^I/ = k 2 (? m ^] for comparison 
the time derivative terms in equation (16) are ~ 7i 2 ^> . Thus, for wavelengths much longer 
than the comoving Jeans length Xj defined by tt/Xj = 7i/c w , the sound wave term may be 
neglected. In the standard cosmology, the Jeans length at z < 100 is less than about 20 
Mpc. 

The entropy source term in equation (16) can also be obtained by comparing separate 
Robert son- Walker spacetimes with slightly different entropies, allowing us to reduce the 
long- wavelength evolution entirely to quadratures, as we show next. 
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2.3. Entropy Perturbations 

Consider two Robert son- Walker spacetimes with identical spatial curvature K but with 
entropies S and S + SS respectively, where SS is a small constant. At a given expansion 
factor a (a 3 plays the role of volume), the density p(a, S) and pressure p(p, S) will differ 
slightly in the two spacetimes. Writing the pressure as p(p(a, S), S), from equations (9) and 
(13) we obtain 

dp\ ( dp\ 
dSJ„ = \dSj 



dp 

d In a 



-3cl(p + p) , — 



+ C -h- 



dp 

dS 



;i7) 



Using equations (9) and (17), we find 

dp\ fp + P 



d 



d In a 



P + P 



Op 
OS 



Using this result, we define a fractional entropy perturbation variable 



a(a, S) 



Sp - c 2 Jp = J^ / dp_ \ 
p + p P + P \dS J 



5S d 



3 d\na [p + p 



dp 
dS 



(19) 



This is a formal result because for a multicomponent imperfect fluid one would not evaluate 
p(a, S) but would instead characterize density and pressure perturbations for the individual 
components, as we shown below in §3.1. For now, we assume that a can be determined and 
we use it to derive a quadrature for the isocurvature modes. 

The scale factor in the second Robert son- Walker spacetime is a(r, K, S+5S) = a(r, K, S)+ 
(da/dS)5S. Changing time variable r — > r + a(r) (with no change in spatial coordinates), 
the line element for the second spacetime becomes 



ds 2 = a 2 (r, K, S) (l + 2Ha + 25S^^\ [-(1 + 2a)dr 2 + d X 2 + r 2 ( X , K)dQ 2 ] . 



(20) 



This line element describes a perfectly homogeneous and isotropic Robertson- Walker space- 
time with entropy S + 5S. However, for appropriate choice of a (i.e., the appropriate coor- 
dinate transformation), it takes precisely the same form as a perturbed Robertson- Walker 
spacetime with background entropy S, given by equation (5), provided that the following 
two conditions hold: 



2^ 



n 



* 



-5S 



d In a 

dS 



Ha . 



(21) 



So far we have assumed nothing about gravity. To proceed further we assume that 
equation (7) is valid, giving 

d\na\ (dT/dS) a ^ K 



OS 



(dT/d\na) K „ 



WS^t(I)/- <*> 
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Combining equations (10), (19), (21) and (22) gives 
1 d fa 2 ^\ fl 1 fO\na 

— OD- 



o? dr V H J dr 



H\ dS 



T,K 



47rGa 2 SS f dp 



3H 2 V OS 



1_ 

n 2 



[ a(a(r', K, S), S)H(t', K, S) dr' . (23) 

This result agrees exactly with equation (16) for long-wavelength entropy perturbations in 
GR with vanishing shear stress potential n. Therefore, in GR the dynamical evolution 
of long-wavelength entropy perturbations in a Robertson- Walker spacetime follows directly 
from the Friedmann and energy conservation equations without requiring the perturbed 
Einstein field equations. Equation (23) implies equation (12) when the entropy perturbation 
vanishes. 



2.4. General Solution for Long- Wavelength Perturbations in GR 

Equations (12) and (23) have a simple exact solution. The homogeneous solution with 
o = (curvature perturbations) is (suppressing the dependencies on K and S, which are 
held fixed) 

h r ^Wjr') n 

HHt>) ' { ' a 2 



*(t)=kMt)+C*.(t), V + (r) = i \ J 2 ,\ ' dr' , <Mt) = 4. (24) 



Here k is the curvature perturbation of equation (11) and it gives the amplitude of the 
"growing mode" of density perturbations; because the lower limit of integration for ^+(t) is 
unspecified one may add any constant multiple C of the "decaying mode" solution \&_(t) = 
7ia~ 2 . The decaying mode is a gauge mode which can be eliminated using the coordinate 
transformation r — ► r+Ca~ 2 . The particular solution with \l/ = , ]/ = 0atr = (isocurvature 
perturbations) and c^O (again suppressing the dependencies on K and S) is 

*(r)=/ [^ + (r)*_(T / )-*+(T / )*_(r)](T(T / )o 2 (T / )dr / . (25) 

Jo 

The general solution of equation (16) with n = A^ = is given by adding equations 
(24) and (25). We see that it is equivalent to integration of 

„. 4 + 1 ,,.-*, + .(i- 2 |*±)-«^ , 26 , 

when the Friedmann equation governs the background expansion. By comparing the evolu- 
tion of separate Robertson- Walker universes we have solved equations (26) by the quadra- 
tures of equations (24) and (25). The solution requires only that the background evolution 
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obeys the Friedmann and energy conservation equations, and that spatial gradient terms 
are negligible. Later we will drop the assumption of the Friedmann equation to obtain 
quadratures for any theory of gravity. 

Although we have set k 2 = —A = 0, the solution is valid for all wavelengths much 
greater than the Jeans length. One simply allows k and SS to depend on wavevector k as 
determined by initial conditions. At short wavelengths (shorter than the Jeans length, for 
example) k and SS, as obtained using equations (11) and (19), depend on time. 



3. Density and Velocity Perturbations and Einstein Constraints 

In addition to the metric perturbations, the density and velocity perturbations of the 
matter can also be reduced to quadratures in the long-wavelength limit assuming local 
energy-momentum conservation. The density and velocity perturbations then obey initial- 
value constraint equations which relate them to the metric perturbations. 

Under the coordinate transformations given in the preceding sections, the value of T° 
does not change to first order in the perturbations. However, comparing the density field 
of two different Robertson- Walker spacetimes at the same coordinate values (r, x l ) gives a 
density perturbation: 4 



S(t) 



P + P 



, d In a \ . „ ( d In a 

Ha + 2n[ — - + SS 



dlnKJ „ V dS , , 



SS (dp 

+ JT^{dsj a 



= 3 ^_ K) + Jl__ f9£\ = 3 (^r _ «) _ 3 [ T <j{t')H{t') dr' + A, (27) 
p + p\dSJ a J 

where A is an integration constant. Equation (27) also follows from energy conservation in 
a perturbed Robertson- Walker spacetime, which gives the perturbed continuity equation 

6 + 3Ha = 3^ + Au (28) 

where u is the velocity potential. For long wavelength perturbations of wavenumber k — > 0, 
A-u = — k 2 u is generally negligible compared with the other terms in the equation. 

Determining the velocity perturbations requires considering spatial variations of the 
velocity potential, which do not exist in a perfectly homogenous Robertson- Walker spacetime. 



4 Note the unusual definition of 8; for a gas of particles it gives the relative number density perturbation 
rather than the relative energy density perturbation. The equations simplify with this choice. 
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Thus we consider a perturbed Robert son- Walker spacetime with metric (14) and energy- 
momentum tensor (15). Evaluating V„T^ = gives the first-order result 



u + (1 - 3c 2 w )Hu = c 2 J + a + <S> + (A + 3K)tt . 
If the shear stress can be neglected, we can integrate to obtain 



u(r) 



[cI(t')S(t') + a(r') + $(r')] W) +p{r')]a\r') dr> ■ 



B 



(p + p)a 4 



(29) 



, (30) 



(P + P)a 4 J 

where B is an integration constant. 

Equations (27)-(30) do not assume the Einstein field equations - they hold for any 
theory of gravity provided consistent with local energy-momentum conservation. However, 
the integral forms of these equations are rarely used. Instead, given the potential \& on long 
wavelengths, the usual procedure in general relativity is to evaluate the density and velocity 
potential using initial-value constraints, of which the Poisson equation is one. That these 
constraint equations are more general will be shown next. 

We define the following combinations of metric and energy-momentum perturbations: 



c 3 



(A + 3K)V-~f(8 + 3Hu) , 
if + 7i$ - 7« , 

* - $ - 377T . 



(31a) 
(31b) 
(31c) 



Here, 5 + 37iu is the number density perturbation on a hypersurface in which T° i vanishes, 
i.e. number the density perturbation in the local fluid rest frame. 

In general relativity each of the constraints vanishes {C\ = C 2 = C3 = 0), however 
we will not assume this to be automatically true. Differentiating the constraints and using 
equations (10), (28) and (29) gives 



-^-(oCi) - (A + 3K - 3 7 )C 2 - H(A + 3K)C 3 = 3 1 (H 2 
a or 



K-H-i)u, (32a) 



2A ( a -c 

a or \7 

7 d 
~H~th 



c 2 C + 7 d ( H2 C 

CwCl + n^\ 7 3 



KC 



n 2 d 



7a 2 dr \H 



,y 



- c 2 w Am - 7 



3 d ,„ ,, s 
° + HTr {H ^ 



A 



7T 



+ (H 2 + K - Ti - 7 )($ - 2*) - [2 7 - 3{r] 2 + K)(l + w)]m . 



(32b) 



Here we set 7 = 47rGa 2 (p + p) and have assumed nothing about gravity. If we assume that 
the Einstein field equations hold, the right-hand sides of equations (32) vanish. If we make 
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the weaker assumption that the background is governed by the Friedmann equation and that 
A\l/ and it (shear stress) can be neglected on large scales so that equation (23) holds, the 
right-hand sides still vanish. In keeping with the previous treatment we also assume \1/ = $ 
so that C3 = 0. With these assumptions we can integrate equations (32) to obtain 

d = mC 2 - A 1 , C 2 = --l T 4(r')7(rV(r') dr' - ^^- , (33) 

fl Jo a, 

where A and B are integration constants. Comparison with equations (27), (30), and (31) 
shows that these are exactly the same initial-value constants obtained from integrating the 
equations of energy-momentum conservation. In other words, they correspond to changing 
5 and u without changing \1/. The third initial- value constant C appearing in equation (24) 
contributes nothing new to the constraints because the pure decaying mode has 5 = 3^ = 
—37iu; the decaying mode can be eliminated by a change of the time coordinate. 

In general relativity, the Einstein field equations give C\ = C 2 = C% = implying 
A — B — 0. However, in other theories of gravity the constraints may be nonzero without 
contradicting equations (32). The initial- value constants follow from 

rp r o 

A = lim — -, B = \ima 4 (p + p)u . (34) 

a^oa 3 (p + p) a->o yi r, v , 

While nonzero values A and B are testable, in principle, by comparing the density and 
velocity fields of galaxies with the gravitational potential implied by gravitational lensing 
or hydrostatic equilibrium in clusters of galaxies, much more stringent tests obtain at high 
redshift, where the A and B terms contribute relatively more to the energy- momentum 
tensor. It would be interesting to place limits on these constants using measurements of the 
cosmic microwave background anisotropy. 

Assuming that A and B are small compared with 5 and a A (p + p)u, respectively, at low 
redshift, and assuming furthermore that the Friedmann equation and energy- momentum 
conservation are valid, then the density and velocity fields in a perturbed Robertson- Walker 
spacetime whose metric takes the form of equation (5) must on large scales obey the same 
constraint equations as in general relativity, namely C\ = C 2 = C3 = 0. 

Although the linearized Einstein field equations give three constraints, the first two 
(Ci = Ci = 0) are initial-value constraints with no dynamical content. If C\ — C 2 — 
on an initial timeslice, energy-momentum conservation combined with the other Einstein 
equations is enough to ensure C\ = C 2 = for all times. The third constraint, C3 = 0, is a 
true dynamical constraint because its time derivative is not forced to vanish as a result of 
the Einstein field equations and energy-momentum conservation. 
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We will see later that initial-value constraints exist not only in GR but in any gravity 
theory yielding a long-wavelength perturbed Robertson- Walker solution with local energy- 
momentum conservation. The key distinguishing features of general relativity are then seen 
to be the Friedmann equation plus the dynamical constraint C3 = 0. 



3.1. Multicomponent Fluids 

Equations (25), (27) and (30) are true quadratures only if <x(r) is known. While equation 
(19) gives a from p(p, S) or p(a, S), the entropy depends on the internal degrees of freedom 
of a fluid. 

Consider a multicomponent imperfect fluid whose density and pressure are described 
by a set of parameters {X} which can vary with position at fixed expansion factor a. We 
replace the single entropy S by as many parameters as are necessary to characterize spatial 
variations in the equation of state of the fluid. For example, a system of non-interacting 
fluids has 



p(a, X) = ^ Pi r i( a ) > p(a,X) = ^WiPir^a) , r;(a) = exp 3/ (l + w^dlna 

L J a 



(35) 



where the pi are independent of a. The Wi generally are not spatially varying, for example 
w = for cold dark matter ("dust") and w = | for a relativistic ideal gas ("radiation"). In 
this case the parameters are the abundances of each fluid component, {X} = {pi,p2, ■ • •}• 
The same procedure will work regardless of what parameters characterize the multicompo- 
nent fluid but the set {X} includes only those that can vary with position at a fixed value 
of a (a playing the role of volume). 

The net equation of state and sound speed parameters are w = p(a,X)/p(a,X) and 
c^ = (dp/ da) /(dp/ da). The entropy perturbation follows similarly to equation (19): 

^ 5X ( dp 2 dp \ 1 de _ ^^ 5X f dp \ 

° = 22—p{dX- C -dx) = -3dhra> < a > X ^ = i:—p{dx) a - (36) 

The symbol e, equalling the number density perturbation at fixed expansion factor a, is 
introduced here to avoid confusion with the number density perturbation 5 measured at 
fixed r. For the example of equation (35), 

° = ^2fi(wi- c 2 w )ti = Y^ fm(ti ~ e) , (37) 

i i 

where 

fe+^iM, f , s * £ = v>,. (38) 

P + V Pi+Pi ^ 
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Here pi = Wipi and e« are both independent of time and fi is the enthalpy fraction for 
component i, normalized so that Yliifi = 1- We see that the entropy perturbation of a 
superposition of ideal gases arises entirely from differences in the particle number density 
of the different species measured at fixed expansion a. For example, for cold dark matter 
(subscript m) and radiation (subscript r) with y(a) = p m /p r , fm = 3y/(3y + 4), f r = 
4/(3?/ + 4), and a = \f r fm(e r — e m ). For a superposition of N fluids there are N — 1 
independent entropy modes corresponding to the N—l independent differences of the number 
density perturbations tj. 

In the general case where the individual components have a time-varying Wi and may 
also have entropy perturbations, 



*=>>[* + {c\ - 4K] =\U [a + c\(e % - e)] , (39) 



where cf = [d(wiri)/da]/(dri/da) and a t are the sound speed squared and entropy pertur- 
bation for each component, respectively. As an example, consider a tightly-coupled plasma 
of photons and nonrelativistic baryons with w = [3(1 + yb)}~ 1 where yb{a) = Pb/Pr (subscript 
b denotes baryons with r for the photons). In this case, w = / r /(4 — f r ) and (? w = |/ r 
compared with w = / r /(4 — f r ) where f r = 4/(3yb + 4). Although this plasma behaves like 
a single fluid, it can have a nonzero entropy perturbation arising from initial variations in 
the baryon to photon ratio, a = \f r fb{^r ~ Q>) where e r — 65 = 5m(n r /n&), the fractional 
perturbation in the photon to baryon number density ratio, is constant in time. Thus, by 
examining the composition of each fluid component, we can write the time-dependent total 
entropy perturbation a in terms of a set of time-independent constants Q. By doing so, we 
are able to fully reduce to quadratures the metric perturbation as shown in §2.4. 

Using the same methods we obtain exact results for the density of individual fluid 
components. Here one must be careful to evaluate the perturbations at fixed (r, x) after the 
coordinate transformations r — > r + a(r) and x ~ * x(l — K ) discussed in §2.1. The density 
perturbation of component % at fixed r is 

Pl (a(r + a,K + 5K,X + 5X),X + SX) - Pl (a, X) 

Oi{T) = ; = 3(W - K) + 6i . (40) 

(P + Ph 

The €i comes from the initial number density perturbation at fixed a while 3(k — \P) is the 
fractional change in volume introduced by shifting the initial hypersurface of constant r. 
Averaging over all components with weights fi gives the net density perturbation 

S( r ) = 3(W - «) + e(r) , (41) 

where e(r) = e(a(r,K,X),X) is given by equation (36) and we are holding K and X fixed 
to lowest order in perturbation theory. By comparing equations (27) and (41) and using 
de/d In a = —3a we see that the constant A has been determined. 
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The velocity potential of individual components follows from equation (30) with sub- 
script % applying to all quantities except a, r, and $. For N uncoupled fluids there are N 
constants Bi. The net velocity potential is the average over the fluids, u = £V fiU u . 

We have succeeded in reducing the long-wavelength density, velocity, and entropy pertur- 
bations of all fluid components to quadratures. This presentation works for both general rel- 
ativity and alternative theories of gravity. In the former case, equations (24) and (25) reduce 
the problem entirely to quadratures specified by a set of constants (k, e±, . . . ejy, -B 1; . . . , B N ). 
We explore next the reduction of the metric perturbation to quadratures for alternative 
theories of gravity. 



4. Quadratures for Alternative Theories of Gravity 

Much of the treatment given so far assumes the validity of the Friedmann equation 
or GR. It is straightforward to redo the calculations of equations (4)— (11) and (22)-(23) 
without making these assumptions. The expansion rate TL can depend only on the scale 
factor a and curvature K of the background Robertson- Walker spacetime and the parameters 
{X} describing the equation of state as discussed in §3.1: 7i = Tt(a,K,X). In this case 
equation (7) is replaced by some function r(a,K,X) = J Q a [a7i(a, K, X)]" 1 da. In addition, 
for an arbitrary theory of gravity the scalar mode has two distinct gravitational potentials, 
equation (14). Repeating the derivations of §§2.1 and 2.3 with the metric of equation (14) 
with $ 7^ ^ and with an arbitrary background expansion rate 7i(a, K, X) gives the following 
result for long-wavelength metric perturbations with 7r = 0: 



1 d fa 2 ^\ _ k d f a\ „, d 

■($-#) = ——( — ) +H- 



o? dr \ H J adr \Ti) d\na 



, di 
2k 






9hlK 'a,X x \—/a,K 



- ^y- 2K fcL"?^(i^L- (42) 

The quantities k and 5X are (for long wavelengths) independent of time. 

Equation (42) provides one differential equation for two functions $ and \1/. Without a 
second relation, provided by a law of gravity, we can determine neither $ nor \1/. 

Even without a full theory of gravity, however, we can still obtain interesting and useful 
results if we assume energy-momentum conservation. In this case we have the quadratures 
(30) and (41). By combining these equations with equation (42) one can derive the following 
results which generalize the energy-momentum constraints C\ and C<i of general relativity 
in the limit of vanishing shear stress (n = 0) and long wavelengths (the spatial Laplacian 
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A = when applied to all perturbations): 

3HB' 



5 + 3Hu 


a\p + p) ' 


(43a) 


y + H$ 


d fn\ H d (HB'\ 
dr \a J a 2 (p + p) da \ a ) 


(43b) 


where 






B\a) = [ a 2 (p+p) 
Jo 


-2k ( 


,d\nK) aX L " V dX ) a , K + 3adr \h)k* _ 


ada 

n +B - 



(44) 
The constant B is essentially the same integration constant as in equation (30). Without 
loss of generality we may set B = in equation (44). 

Equations (41)-(43) are not all independent. Any three of the four imply the remaining 
one. Thus, we may describe the perturbations by equations (41) and (43). To reduce the 
system fully to quadratures requires one more equation equivalent to a relation between $ 
and ^. 

It is worth emphasizing that the reduction to quadratures is based on only a few as- 
sumptions: 

1. The spacetime is nearly Robertson- Walker with perturbation amplitude small enough 
for linear perturbation theory to apply. 

2. Local energy-momentum conservation holds: V^T^ = 0. 

3. In the conformal Newtonian gauge, spatial gradients of fields are small enough to be 
neglected in all equations of motion. 

4. Shear stress perturbations are negligible. 

5. A theory of gravity provides some relation between $ and ^ and possibly the matter 
perturbations. This assumption is needed only for a complete reduction to quadratures. 



There are no restrictions on the equation of state of matter and radiation fields. There are 
no restrictions on the geometry of perturbations; in particular, no assumption of spheri- 
cal symmetry. The neglect of spatial gradients becomes invalid on scales over which non- 
gravitational forces act so as to modify local energy-momentum conservation. By definition, 
spatial gradients of pressure are important and modify our results on scales less than the 
Jeans length. 
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The constraints (43) are particularly useful when B' = 0. For any theory having a flat 
Robert son- Walker solution (including GR), lim^^ ( , 9hi7Y/(91ni^) = 0. When there are no 
entropy perturbations, e = 0. Thus, in a flat universe with initially isentropic perturbations, 
the only possible contribution to B' comes from the J2x SX (dlnH/ din X) term. If the only 
dependence of Ti on equation of state parameters X is through the density p, then 



x v 



dlnH 
dX 



e(p + p) 



a,K 



dlnH 
dp 



(45) 



which vanishes for isentropic initial conditions e = 0. Thus, for a broad class of theories, 
B' — for curvature perturbations in a flat universe implying the long- wavelength constraints 
5 + 3?iu = 0, \1/ + H$ = 7« where 7 = —ad(7i/a)/dr for K = 0. From equation (41) we also 
have 5 = 3(^/ — k) for long wavelengths and isentropic initial conditions. In this case we have 
three relations between the four time- dependent functions ($, \l/, S, u). In the general case we 
have an additional function, B', whose specification requires a theory for the homogeneous 
expansion. 

Thus, without knowing anything about the underlying gravity theory except that it 
is consistent with local energy-momentum conservation and that it admits a Robertson- 
Walker solution, we have deduced initial-value constraints similar to those of GR under the 
conditions of long wavelengths and negligible shear stress. One is not free to specify arbitrary 
metric and matter perturbations on long wavelengths, at least under the assumption that 
long-wavelength perturbations evolve like separate universes. 

However, the initial-value constraints (43) combined with equations (36) and (41) are 
insufficient to fully determine the growth of perturbations. We need a third constraint, a 
relation between $ and \l/ given by some theory of gravity. Such a dynamical constraint would 
allow equation (42) to be integrated and, with the initial- value constraints, reduce ($, \1/, 5, u) 
fully to quadratures for long-wavelength perturbations with vanishing shear stress. 

The simplest case of a dynamical constraint on the gravitational fields is $ = ^ for 
71 = as in GR (although this is more general than GR since we do not assume the Fried- 
mann equation). Under this assumption, equation (42) may be integrated completely giving 
(suppressing the dependencies on K and X and using a as the time variable) 



#(a) =K* + (a)+J](0:)*x(a)+C#_(a) if $ = ^ 



(46) 



x 



where 



*_ 



n 

7^2 



1- 



dlnH 
d In a 



-2 



K,X 



dlnH 
d\nK 



a,X 



ada 



VJa) 



Tia 



(47) 
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The potential has been reduced entirely to quadratures depending only on the expansion 
rate 7i as a function of expansion, curvature, and equation of state parameters. 

Considering again a flat universe with isentropic initial conditions and assuming equation 
(45) holds, but now with the added condition $ = \I/, we see that ^ = k^> + is determined 
completely by the expansion rate 7i(a,0,X) with constant values of the equation of state 
parameters X. The density and velocity potential perturbation follow from S = 3(\l/ — k) — 
—37iu. In this case the growth of long- wavelength perturbations is completely determined 
by the expansion history 7i(a) with fixed curvature (K = 0) and composition (5X = 0). 

For example, in the DGP brane- world theory of Dvali et al. (2000), the Friedmann 
equation is modified by replacing p by (VP + Pr c + \fp~rl) 2 with p Tc being a constant. If 
$ = \I/, equation (11) remains valid for the long- wavelength curvature mode provided that 
7 changes to 

i = H 2 + K - H = I (1 + w) (H 2 + K) (^l£±PjJL \ . (49) 

2 \VP + Prc + VKJ 

The isocurvature mode is also modified slightly, with 7 — ► 7' p/(p + p Tc ) in the last line of 
equation (23). Equations (27) and (30) are unchanged. 

The solutions implied here are only valid as k — ► so that spatial gradient terms 
may be neglected in the equations of motion. The DGP model has a length scale r c = 
(327rGp rc /3) -1 / 2 comparable to -f^o" 1 . The long-wavelength limit becomes kr c <^i 1 suggesting 
that the formal solution given here is not applicable for wavelengths shorter than the present 
Hubble length. Other theories of gravity will give different ranges of applicability. In general, 
any long-range departures from general relativity will modify the evolution of patches of a 
Robertson- Walker spacetime. 

The results of this section highlight the importance of testing the dynamical constraint 
C3 = ty — $ — 377r = in general relativity. The relation between $ and \1/, combined with 
the background expansion rate 7i(a, K, X), is the key to long-wavelength perturbations of a 
Robertson- Walker cosmology. 
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5. Comparing Geometric and Dynamic Methods for Testing Gravity and 

Dark Energy 

As we have seen, the evolution of long- wavelength perturbations requires two ingredients: 
(1) a relation between the two gravitational potentials $ and \P (possibly involving auxiliary 
fields), and (2) a generalization of the Friedmann equation for the background expansion 
rate. Without the first ingredient, the long-wavelength perturbations cannot be predicted. 
Let us assume that a dynamical constraint exists relating $ and \l/ so that equation (42) can 
be integrated. Then the evolution of perturbations is determined by the expansion rate of 
the universe, the very same information that determines the redshift- distance relation used 
by geometric methods. Does this mean that dynamic and geometric methods measure the 
same thing? 

Not necessarily. First, the dependence of perturbation (dynamic) methods on the expan- 
sion rate necessarily depends on the value of $ — \1/ which may be nonzero for theories of grav- 
ity other than GR. Second, dynamic methods depend not only on H(a,K,X) and its vari- 
ation with a but also on (d\nH/d\nK) aX (for curvature perturbations) or (<91n7Y/<9X) a x 
(for entropy perturbations; here X are parameters describing the fluid composition or equa- 
tion of state). Geometric methods, by comparison, depend only on 7i(a,K,X) with fixed 
K and X. If the universe is nonflat or there exist entropy perturbations, dynamic methods 
have the potential to reveal information about the dependence of 7i on K and X which 
is not present in the expansion history for one universe with fixed K and X . Finally, the 
quadrature results for perturbations are valid only in the limit of long wavelengths; long- 
range forces different from Einstein gravity may lead to spatial gradient terms that modify 
the quadratures on presently observable scales. 

Geometric tests are generally cast in terms of the luminosity distance or angular- 
diameter distance but in fact depend strictly on the redshift-distance relation for radial 

null geodesies, 

r dz' 

X (z) = r -r(a,K,X) = J ^, (50) 

where To = r(l, K,X) and a = (1 + z)^ 1 . The redshift-dependence of the long-wavelength 
curvature perturbation \P+(z) requires specifying a theory of gravity. If $ = ^/, for both 
general relativity and DGP brane worlds the same result holds on very long wavelengths: 5 



*_ 



[l + z)H(z) 



K 



+ 



1 



dH 



H 2 (z') ' (1 + z')H(z') dz' 



dz' 
H(z' 



(51) 



5 This does not mean that the two theories predict identical perturbations on scales larger than the Jeans 
length. The DGP model has a much longer length scale, r c , below which equation (51) may be invalid. 
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More generally, let us make no assumptions about $ — \I/ but consider a flat universe 
with negligible shear stress and initially isentropic curvature fluctuations as predicted by the 
simplest inflationary universe models. In this case, equation (42) reduces to the following 
relation between H(z), ^{z), and &{z) for long wavelengths: 

$(*) _<** , (K-*) dH 

1 + z dz H dz 

Here k is independent of redshift on large scales although k, \I/, and $ will vary with position 
or spatial wavenumber. In principle, measurements of H(z) from geometric methods and 
ty(z) from perturbations could (for isentropic perturbations of a flat universe with negligible 
shear stress) determine $(z) up to an additive term proportional to dlnH/dln(l + z) (since 
the relation between k and \1/ is unknown if the theory of gravity is unspecified). In particular, 
such measurements could enable GR or alternative theories of gravity to be tested without 
any assumptions about the density and pressure of mass-energy in the universe. 

This kind of test is difficult to imagine carrying out because of the difficulty of measuring 
^(z). 6 Alternatively, if one assumes that GR is valid, dynamic methods can be used to 
provide additional constraints on dark energy because the geometric and dynamic methods 
have a different dependence on the equation of state of dark energy. 

Geometric measurements of x( z ) a ^ a given redshift depend on the expansion history 
at smaller redshifts; dynamic measurements of ^(z) depend on H(z) at higher redshifts. 
Because the effects of dark energy typically decline with increasing redshift, perturbation 
methods are at a disadvantage at high redshift but may be superior to geometric methods 
at low redshift. 

To assess their relative merits for measuring dark energy (assuming GR is the correct 
gravity theory) the methods were compared using a flat model with Q m = 0.3 and a dark 
energy component with constant w independent of z. While unphysical, this model is com- 
monly used to compare theories with data. The logarithmic derivatives of equations (50) and 
(51) with respect to w at fixed z were evaluated numerically to produce the results shown 
in Figure 1. A large value of \dlnx/dw\ indicates that the geometric method is relatively 
powerful — a given measurement error in In x translates into a smaller uncertainty in w for 
larger values of \dlnx/dw\. The same is true for dynamic methods using ln , I' + . 

As expected, perturbation methods are insensitive to the dark energy parameter w at 
high redshift. At low redshift, dark energy dominates more rapidly for smaller w leading to 



Observations of density perturbations, e.g. from galaxy redshift surveys, are easier than observations of 
the gravitational potential perturbations. The relation between density perturbations and ^f depends on the 
theory of gravity and may require specifying more than just H(z). 
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Fig. 1. — Sensitivity of geometric methods (left panel, characterized by the comoving dis- 
tance x to redshift z) and dynamic methods (right panel, characterized by the curvature 
perturbation \l/ + at redshift z) to the equation of state parameter w. A flat model with 
Q m = 0.3 was assumed. Curves are shown for w = —0.6, —0.7, . . . , —1.1. Dynamic (pertur- 
bation) methods are insensitive to dark energy at high redshift but are more sensitive than 
geometric methods at low redshift. 

a greater suppression of linear growth hence a larger variation with w. Geometric methods, 
on the other hand, are insensitive to w at small redshift where the lookback time is much 
less than the Hubble time. For the standard Friedmann equation and w ~ — 1, dynamic 
methods are more sensitive than geometric methods for z < 0.2. At high redshift the 
comoving distance remains sensitive to w despite the declining importance of dark energy 
because x( z ) i s an integral over the past lightcone. The expansion history at low redshift 
affects the redshift-distance relation at high redshift. 

The theoretical sensitivities shown in Figure 1 must be combined with the measurement 
uncertainties of the two methods before a reliable estimate can be made of their relative 
merits. The perturbation methods face a severe challenge - - to achieve a discriminating 
power of 0.1 in w requires 2% accuracy in measurement of \I/+ at low redshift. 



6. Summary and Outlook for Cosmological Tests of Gravity 

Long- wavelength cosmological perturbations involve at least two metric perturbations 
($, \l/) and two matter perturbations (S,u) (density and velocity potential, all in conformal 
Newtonian gauge). One might think that, without a theory of gravity, no predictions can be 
made about relations between these variables. 
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In fact, we have shown that, with minimal assumptions, on large scales there are auto- 
matically three independent relations between ($, \l/, 5, u). For an arbitrary theory of gravity 
these are any three of the four equations (41), (42) and (43a,b). Equations (43) enforce lo- 
cal energy-momentum conservation. They generalize the initial- value constraints of general 
relativity. GR also has a dynamical constraint on $ — \I/ enabling a complete reduction 
to quadratures. These relations depend on the expansion rate of the background space- 
time and on the dependence of the expansion rate on spatial curvature (if K ^ 0) and on 
the composition or entropy of the matter filling space (if there are composition or entropy 
perturbations). 

Thus, the key ingredients needed to specify long- wavelength perturbations are (1) a 
relation between the two gravitational potentials $ and \1/, and (2) a relation between the 
expansion rate, density, and pressure (e.g. the Friedmann equation). Given these two ingre- 
dients, we have shown how to reduce ($, ^,5, u) to quadratures by introducing conserved 
curvature and entropy variables. Explicit expressions for the time-dependence of the metric 
perturbations were given for Einstein gravity, equations (24) and (25) or (46)-(48). 

In some circumstances (e.g. isentropic fluctuations in a flat universe with $ = \l/) 
the long- wavelength growth of perturbations is determined completely by the expansion 
history H(z). In such cases measurements of perturbation growth cannot be combined 
with geometric measurements (e.g. supernova distances) to discern whether dark energy is 
a new form of mass-energy or a failure of the Friedmann equation. A loophole exists in 
this argument if $ 7^ ^ or there are significant long-range non-gravitational forces. For a 
flat universe with curvature fluctuations, the argument can be inverted to yield information 
about *&(z) from measurements of H(z) and ^(z) thereby providing a test of gravity theories 
independently of assumptions about the energy density and pressure. 

The long-wavelength limit corresponds to wavelengths larger than any relevant spatial 
scales so that spatial gradients may be neglected in the equations of motion. Pressure forces 
modify the dynamics on scales smaller than the Jeans length; the dark energy may plausibly 
be a scalar field with a Jeans length comparable to the Hubble length. By measuring the 
wavelength-dependence of the linear growth rate on scales greater than 1 Gpc one might 
measure the time-dependence of the dark energy Jeans length and thereby constrain its 
intrinsic properties. This measurement is exceedingly difficult because the amplitude of 
density perturbations in the dark energy is expected to be less than about 10 -4 . 

In addition to providing tests of general relativity, perturbation measurements can pro- 
vide constraints on dark energy under the assumption that general relativity is valid. By 
comparing the dependence of curvature perturbations and the redshift- distance relation on 
the dark energy equation of state parameter we verified quantitatively the expected con- 
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clusion that perturbation methods are most useful at small redshift when the accelerated 
expansion begins to suppress curvature perturbations. A changing gravitational potential 
generates cosmic microwave background anisotropy through the integrated Sachs- Wolfe ef- 
fect. Measurement of this effect (Padmanabhan et al. 2005) has the potential to further 
constrain dark energy (Pogosian et al. 2005). 

The analysis in this paper shows the importance of testing the equality of the two 
Newtonian gauge gravitational potentials, $ = \1/ in equation (14). While this equality 
holds in general relativity (in the absence of large shear stress), it may not be true for other 
theories of gravity. In addition to tests combining geometric and perturbation methods using 
equation (52), one can in principle measure $ — ^ by comparing the deflection of light by 
gravitational lenses (an effect proportional to ^ + $) with the non-relativistic motion of 
galaxies (an effect proportional to $). A similar test exists on solar-system scales (or in 
binary pulsars) where the deflection (or Shapiro delay) of light is compared with Newtonian 
dynamics. Thus, for testing modified gravity as an alternative to GR it is important to 
extend tests of the post-Newtonian parameter 'Jppn = ^/$ (Will 2006). Stringent limits on 
\1ppn — 1| apply on solar system scales. It would be worthwhile to improve the limits on 
Mpc and larger scales by combining weak gravitational lensing and galaxy peculiar velocity 
measurements or by adding 7ppjv to the parameters used in analyzing cosmic microwave 
background anisotropy. 

I thank Eric Linder for helpful comments. This work was supported by the Kavli 
Foundation and by NSF grant AST-0407050. 
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